Magnetic properties and thermal entanglement on a 
triangulated Kagome lattice 



N S Ananikian^, L N Ananikyan^, L A Chakhmakhchyan^'^ and 
A N Kocharian^'^ 

-'^Yerevan Physics Institute, Alikhanian Br. 2, 0036 Yerevan, Armenia 

^ Yerevan State University, A. Manoogian 1, 0025 Yerevan, Armenia 

^ Department of Physics, Cahfornia State University, Los Angeles, CA 90032, USA 

E-maih ananikSyerphi . am 

Abstract. The magnetic and entanglement thermal (equilibrium) properties in spin- 
1/2 Ising-Heisenbcrg model on a triangulated Kagome lattice are analyzed by means of 
variational mean-field like treatment based on Gibbs-Bogoliubov inequality. Because 
of the separable character of Ising-type exchange interactions between the Heisenberg 
trimers the calculation of quantum entanglement in a self-consistent field can be 
performed for each of the trimers individually. The concurrence in terms of three 
qubit isotropic Heisenberg model in effective Ising field is non-zero even in the absence 
of a magnetic field. The magnetic and entanglement properties exhibit common 
(plateau and peak) features observable via (antferromagnetic) coupling constant and 
external magnetic field. The critical temperature for the phase transition and threshold 
temperature for concurrence coincide in the case of antiferromagnetic coupling between 
qubits. The existence of entangled and disentangled phases in saturated and frustrated 
phases is established. 
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1. Introduction 

Geometrically frustrated spin systems exhibit a fascinating new phases of matter, 
a rich variety of unusual ground states and thermal properties as a result of zero 
and finite temperature phase transitions driven by quantum and thermal fluctuations 
respectively P-E] . The quantum statistical studies of perplexing physics due to quantum 
fluctuations, geometric frustration, competing phases and complex connectivity in 
electron and spin lattices still remain a challenging theoretical problem. Some prominent 
properties related to the frustrated models, such as transition from disordered quantum 
spin liquid into spin insulator with broken translational symmetry similar to the Mott- 
Hubbard insulator at half filling with antiferromagnetic plateau behavior, seen in 
magnetization versus magnetic field at low-temperature, have been recently intensively 
studied both experimentally [TUTU] and theoretically [TTHT3] . Furthermore, the (non- 
bipartite) frustrated local geometry in two and three dimensions can provide also 
insight into the magnetism of strongly correlated electrons in small and large bipartite 
and non-bipartite (frustrated) Hubbard lattices also away from half filling [H]. The 
efforts aimed at better understanding of the aforementioned phenomena stimulated 
an intensive search of two-dimensional geometrically frustrated topologies, such as 
recently fabricated metalo-organic compound piperazinium hexachlorodicuprate jl5] . 
AgNi02 triangular magnet [IE], crystal samples of KFe3(OH)6(S04)2 (as an ideal 
Kagome lattice antiferromagnet) [17], alternating Kagome and triangular planar layers 
stacked along the [111] direction of the pyrochlore lattice [18], etc. Besides, studies 
of inorganic molecular materials with paramagnetic metal centers connected in the 
crystal lattice via superexchange pathways, strongly frustrated by their geometric 
arrangement, has also been implemented [19]. From this perspective, the most 
interesting geometrically frustrated topologies are the magnetic materials in form of 
two-dimensional isostructural polymeric coordination compounds Cu9X2(cpa)6 • nH20 
(X = F, CI, Br and cpa=carboxypentonic acid) [20l - [22] . The magnetic lattice of these 
series of compounds consists of copper ions placed at two non-equivalent positions, which 
are shown schematically as open and full circles in figure [H Cu^"*" ions with a square 
pyramidal coordination (a-sites) form equilateral triangles (trimers) which are connected 
one to another by Cu^^ ions (monomers) with an elongated octahedron environment 
(b-sites) forming the sites of Kagome lattice. This magnetic architecture, which can 
be regarded as triangulated (triangles in triangles) Kagome lattice, is currently under 
active theoretical investigation [23|[2^ . 

The spin- 1 Ising model on the triangulated Kagome lattice has been exactly solved 
in [25]. However, in its initial form theory fails to describe the properties of the 
aforementioned compound series, since it entirely neglects quantum fluctuations firmly 
associated with a quantum nature of the paramagnetic Cu^+ ions having the lowest 
possible quantum spin number 1/2. Further extension to the Ising- Heisenb erg model 
by accounting for quantum interactions between Cu^"*" ions in a-sites (with quantum 
spin number 1/2) in the limit when monomeric 6-spins having an exchange of Ising 
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character, provides much more richer physics and displays essential features of the copper 
based coordination compounds [26|l27]. The strong antiferromagnetic coupling has been 
assumed for Jaa between trimeric a sites, while there exists a weaker ferromagnetic 
exchange Jab between the trimer a- and monomer 6-sites at the ratio \JablJaa\ ~ 0,025 

[2H]. 

Entanglement is a generic feature of quantum correlations in systems, which 
cannot be quantified classically [29l|30]. It provides a new perspective to understand 
the quantum phase transitions (QPTs) and collective phenomena in many-body and 
condensed matter physics. This problem, which has been under scrutiny for nearly two 
decades, has attracted much attention recently |3TH3l] . 

A key observation is that quantum entanglement can play an important role in 
proximity to the QPTs, controlled by quantum fluctuations near quantum critical 
points |35] . A new line research points to a connection between the entanglement of a 
many-particle system and the existence of the QPTs and scaling [361438] . Thermal 
entanglement was detected by experimental observations in low dimensional spin 
systems formed in the compounds Na2Cu5Si40i4 |39], CaMn2Sb2 [IQ], pyroborate 
MgMnBsOs the warwickite MgTiOBOg [41j, KNaMSi40io (M=Mn, Fe or Cu) [32] 
and metal carboxylates from magnetic susceptibility measurements [13]. Entanglement 
properties for few spins or electrons can display the local intrinsic features of large 
thermodynamic systems and can be suitable for calculations of basic magnetic quantities 
and various measures in entanglement associated with QPTs. The basic features of 
entanglement in spin-| finite systems are fairly well understood by now (see e.g. [SUlS]), 
while the role of local cluster topology and spin correlations in thermodynamic limit 
still remain unanswered. There are some approximate methods, such as mean-field like 
theories based on the Gibbs-Bogoliubov inequality, that one can invoke to deal with 
the cases like this, aimed at better understanding of different physical aspects [36]. 
This method can also be applied for studying thermal entanglement of many-body 
systems [1?]. In spite of the method not being exact, it is still possible to observe 
regions of criticality [33 • 

In the case of the triangulated Kagome lattice each a-type trimer interacts with its 
neighboring trimer through the Ising-type exchange, i.e. classical interaction, therefore 
the states of two neighboring a-sublattices become separable (disentangled) [29j. Thus 
we can calculate concurrence (a measure of entanglement [3^), which characterizes 
quantum features, for each trimer separately in self-consistent field. The key result of 
the paper is concentrated on the comparison of specific (peak and plateau) features in 
magnetization, susceptibility, specific heat and thermal entanglement properties in the 
above mentioned model using variational mean-field like approximation based on Gibbs- 
Bogoliubov inequality. We will demonstrate how the order- disorder phase transition 
temperature is relevant to the the threshold temperature for vanishing of entanglement. 

The rest of the paper is organized as follows: in section |2] we introduce the 
Ising-Heisenberg model on the triangulated Kagome lattice and provide a solution in 
variational mean-field like approximation based on Gibbs-Bogoliubov inequality. The 
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magnetic properties of the model are investigated in section El The basic principles 
for calculation of concurrence as a measure of entanglement and some of the results on 
intrinsic properties are introduced in section HI In section [5] we present the comparison 
of magnetic properties and thermal entanglement. The concluding remarks are given in 
section |6l 

2. Isotropic Heisenberg model on triangulated Kagome lattice 

2.1. Spin-^ Ising- Heisenberg model 

We consider the spin-| Ising-Heisenberg model on triangulated Kagome lattice (TKL) 
(figured]) consisting of two types of sites (a and b). Since the exchange coupling between 
Cu^^ ions are almost isotropic, the application of the XXX Heisenberg model is more 
appropriate. There is a strong Heisenberg Jaa exchange coupling between trimeric sites 
of a type and weaker Ising-type one (Jab) between trimeric a and monomeric b ones. 
Thus, the Kagome lattice of the Ising spins (monomers) contains inside of each triangle 
unit a smaller triangle of the Heisenberg spins (trimer). The Hamiltonian can be written 
as follows: 

2N 

H = JaaY. S»S» - Jab J^i^Tk " ' H 3[(5^)| + ^{3^% (1) 

(Lj) {k,l) i=l 

where S*^ = {S^, Sy, S"^} is the Heisenberg spin-i operator, is the Ising spin. Jaa > 
corresponds to antiferromagnetic Heisenberg coupling and Jab > to ferromagnetic 
Ising-Heisenberg one. The first two summations run over a — a and a — b nearest 
neighbors respectively and the last sum incorporates the effect of uniform magnetic 
field (we have assumed that the total number of sites is 3A^). 




Figure 1. A cross-section of TKL structure. Solid lines represent the intra- 
trimer Heisenberg interactions Jaa, while the broken ones label monomer-trimer 
Ising interactions Jab- The circle marks k-th cluster (Heisenberg trimer). S^, 
presents the Heisenberg and S^, the Ising spins. 
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2.2. Basic mean-field formalism 

Here we apply the variational mean-field like treatment based on Gibbs-Bogoliubov 
inequality |50] to solve the Hamiltonian ([T]). This implies that the free energy (Helmholtz 
potential) of system is 

F<Fo+{'H- ■Ho)o, (2) 

where H is the real Hamiltonian which describes the system and T/q is the trial one. F 
and Fq are free energies corresponding to T-L and T-Lq respectively and (...)o denotes the 
thermal average over the ensemble defined by "Ho- Following [23j we introduce the trial 
Hamiltonian in the following form: 

"Ho = "Hco, (3 a) 

k&rimers 

3 

1=1 

In this Hamiltonian the stronger quantum Heisenberg antiferromagnetic interactions 
between a-sites are treated exactly, while the weaker Ising-type ones between a- and 
6-sites are replaced by self-consistent (effective) fields of two types: ja and 7^. The 
variational parameters 7^, 7?, and Xaa can be found from the Bogoliubov inequality 
after minimizing the RHS of ([2]). Thus our focus is on the cluster depicted in figured] 
Each of the 6-type spins belongs to two of such clusters simultaneously. Hence only 
the half of each 6-spin belongs to one cluster and therefore there are 3-1/2 b-type 
spins in it. Consequently the total number of (a— and b—) spins in the cluster will be 
3-1-3/2 = 9/2 (there are 3 a-type spins in the triangle). Inequality ([2]) can be rewritten 
for the described cluster: 

f <fo + {'Hc-'Hjo, (4) 

where Tic is the real and 7/co the trial Hamiltonian of the cluster, / and /o free energies 
of the cluster defined by Tic and "Hcq respectively. Using the fact that in terms of ( l3B 
and are statistically independent, one obtains (S" ■ S^)o = (S")o ■ (S'')o. Besides, 
taking into account that {{S^)"')o = ma (single a-site magnetization), {{S^Y)o = nib 
(single 6-site magnetization), we obtain the following expression: 



f < fo + {Jaa - Aaa)(S^^S^2 + ^fe^fca + ^fci^fcs) 







Now, by minimizing the right-hand side of inequality (|5]) with respect to 7^, 7^ and 
and using = -3m., f| = -3/2m,, |£- = (S^^S^^ + S^^S^3 + S^^ 8^3)0, 
we determine the variational parameters in the form: Xaa = Jaa, la = '^Jab^b + H , 
7b = iJabrria + H. Parameters 7. and 7f„ which have a meaning of a magnetic field, 
are interconnected, which is the consequence of its' apparent self-consistency. The 
Hamiltonian Tio was chosen to be exactly solved. One finds that "Hco can be divided 



into two parts corresponding to a- and &-type variables: 



i=l 



(6) 



i=l 



Each of Hamiltonians "H^j^ and {Til^^Y can be solved separately (the variables have been 
separated). The eigenvalues of "H^^ are: 
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^4 = ^ (3A„„ + 27„) ; E, = Ee = ^ (-3A„„ - 27J ; 



(7) 



1 3 

-S7 = - (SXaa — 27a) ; -^'s = ^ (^aa 

and the corresponding eigenvectors given by 



27a) 



1^3) 
1^4) 

1^7) 



|000) 

i=(g|001) + ?>10) + |100)) 



(g2|001) + g|010) + |100)) 
(|001) + |010) + |100)) 



a/3 

1 

V3 

-^(g|110) + g^|101) + |011)) 



(8) 



a/3 



(g2|110) + g|101) + |011)) 



-^(|110) + |101) + |011)) 



iV's) = |111), 

where q — e*^"^/^ (these eigenvectors should be also the eigenstates of cyclic shift operator 
P with eigenvalues 1, q and g^, satisfying q'^ + + 1 = 0). 

The partition function Zq^ of the trimer in mean-field approximation is: 



Zo. = i^^M-Ek/T) = e-^ [cosh (^) 
fe=i L V / 



+ 



(9) 



2e^2T''' cosh 



V2r 



+ cosh (|; 



Consequently the free energy of a-triangle will be: 

3Aaa 



/o„ = -TlnZo„ = 



Tin 



cosh (1^) + 2e^ cosh (^) + (10) 



cosh (|,)_ . 
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Since the CH^^Y describes only half a particle (b-type spin), the Hamiltonian of one b- 
type spin will be 2 ■ (T-Ll^y. Hence, following the technique described above one finds the 
partition function Zq^ and free energy /o^ of a b-type spin in the adopted approximation: 



"Ob 



-Tin 



2 cosh ( , 
\2TJi 



'11a) 



[lib) 



As already mentioned, in terms of the trial Hamiltonian a- and 6-type spins are 
statistically independent [see equations (3) and (E])]. Besides, b— type spins do not 
interact with each other. Therefore the partition function /o of the cluster in mean-field 
approximation reads: 

3 

/o = /oa + 2 ■ fob- (12) 

Consequently the free energy of the cluster fcB in the mean-field approximation based 
on the Gibbs-Bogoliubov inequality will be: 

3A„ 



fcB — fo+ {T^c — '^co)c 



+2 cosh { ^ 
V2T 



4 

2 cosh 



GJabmanif, - T 



In 



4e cosh ( I 
\2TJ 



37a 

2T 



3, 
+-\n 
2 



2-h(^)} 



(13) 



In equation f[T3|) we have used the values of variational parameters 7^, % and Xaa- 
Besides, due to the fact that there are 9/2 spins in the cluster and therefore totally 
2N/3 clusters {Fgb = 2A^/3 ■ /gb), we obtain: 



fGB 

N 



Xn 



+ AJabmairib - 2T 



In < 4e 2T cosh 



\2TJ 



+ 2 cosh 



37a 

2T 



1 



In < 2 cosh 



(-) 
\2TJ 



As for defined above a- and 6-single site magnetizations we obtain: 

1 dfo^ _ 1 3 sinh (1^) + 2e^ sinh (^) + sinh (^) 



rrib 



3 <97a 
<97fe 



6 cosh (1^) + 2e^ cosh (^) + cosh (^) 



- tanh ( 
2 \2T 



(14) 



(15a) 



(156) 



Notwithstanding of simplicity and the fact that the effective (self-consistent) 
field in zero magnetic field {H = 0) overestimates ferromagnetic correlations, it is 
still particularly useful for detection of spontaneous breaking SU(2) symmetry and 
possible temperature driven transitions in the frustrated spin systems (see section \3A}i . 
However, we also find that the strong quantum fluctuations, existing in the isotropic 
Heisenberg model in the absence of ferromagnetic type Ising term at if = can 
restore the broken symmetry by providing stability to disordered spin-1/2 liquid state 
in frustrated geometries. Moreover, in general the presence of magnetic field [H 7^ 0) 
suppresses the spin fluctuations and makes the self-consistent results more reliable and 
accurate [51]. Therefore, the equations (7), (8) and (13)-(15) with magnetic field are 
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quite sufficient for understanding some intrinsic relationships between magnetic and 
entanglement properties [H] that naturally emerge in the Ising-Heisenberg model when 
one is complying with the variational mean-field like procedure. 



3. Magnetic properties 



3.1. Magnetization 

The results of the previous section can be used for investigation of the magnetic 
properties of the model. Here we are interested in the sublattice a properties, which, 
however, depend on parameters describing &-type spins. It is convenient to introduce 
a new (ratio) parameter a = Jab/Jaa- The magnetization curves can be found by 
solving numerically the transcendental equations fllSal) and fll5fe|) . The magnetic field 
dependence of the magnetization per atom is plotted in figure [2] at a = 0.025. We find 
that in the absence of magnetic field, the ordered ferromagnetic phase with spontaneous 
magnetization per site rua is a stable ground state for all \Jab/Jaa\ in spite of the high 
geometric frustration caused by the non-bipartite structure and antiferromagnetic intra- 
trimer interaction. 



(a) 



(b) 



K 0.0 





H 





H 



Figure 2. Single a-site magnetization versus external magnetic field H for 
(a) T = 0.7 and (b) T = 0.01 for a = 0.025, Jaa = 1- The magnetic field and 
temperature are in relative units oi ks = 1. 

At relatively high temperatures the magnetization curve in figure [2](a) shows a 
monotonic behavior versus magnetic field with a full saturation at strong magnetic 
field. Upon decreasing the temperature a new partially saturated phase emerges in 
form of the (spin) plateaus shown in figure [2t^b), which can be associated with staggered 
magnetization or short range antiferromagnetism (AF) in frustrated Kagome geometry. 
Indeed, the appearance of plateau in magnetization curve at rria = 1/6 can be explained 
as stability of trimeric a sites in 'f^l configuration. Thus, at rather low temperatures, 
the magnetization shows the finite leap across a plateau at rria = 1/6 at infinitesimal 
magnetic field and below the critical field for full saturation by fiipping a down spin. 

In figure [3] we also show the temperature dependence of the magnetization in 
equilibrium. As one can see from figure |3]^a), in the absence of the external magnetic field 
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the magnetization tends gradually to zero near the second-order transition temperature 
Tc between ordered (m^ ^ 0) and disordered (m^ = 0) phases. Hence, the magnetization 
rua can be expanded into series near the critical temperature of second-order phase 
transition point: 

■nia = aula + -|- cm^ + ... (16) 

The critical temperature corresponding to the second-order phase transition can be 
found from the condition a = 1, 6 < 0. In particular, for the case Jaa = 1 and a = 0.025, 
Tc = 0.0102062 (in relative units). At zero temperature in the absence of field we find 
unsaturated spontaneous ferromagnetism with = 1/6 as a stable ground state. 





0.015 

T 



Figure 3. Single a-site magnetization nia per atom versus temperature T for 
Jaa = 1, a = 0.025 and (a) H = and (b) different non-zero values of H (the 
inset shows details for the case H = 0.9 at low temperatures). 

The magnetization in equilibrium as a function of temperature T at non-zero 
magnetic field [S2] is plotted in figure El^b). There are two distinct magnetic field regimes 
corresponding to = 1/6 and saturated zero-temperature magnetization, = 1/2. 
While fixed magnetic field H is less than the saturation magnetic field value, we deal with 
TTia = 1/6 regime. If we continue increasing the value of H, at the saturation magnetic 
field the magnetization jumps into nia = 1/2 regime [H = 1.5, 2.0 in figure |3]^b)]. There 
can be also seen a short plateau at nia = 1/6 in the temperature dependence [H = 0.9 
in figure |3]^b)]. The peaks in the case of low magnetic fields arise due to the frustration 
effects. 



3.2. Susceptibility 

We define the magnetic susceptibility Xa as 



Xa 



drrin 



OH ■ 

First we examine the zero-field susceptibility, which is introduced as follows: 



drrir, 



ao 



OH 



(17) 



:i8) 



The temperature dependence of Xm in equilibrium is plotted in figure Si 
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Figure 4. Zero-field susceptibility Xag versus temperature T for Jaa = 1 and 
a = 0.025. 

The zero-field susceptibility Xao diverges at the critical temperature Tc which is a 
signature of the second order phase transition discussed earlier (see section [XT]) . The 
temperature dependence of susceptibility Xa at if 7^ [53] is presented in figure [5]^a). 



(a) (b) 




Figure 5. Susceptibility Xa versus (a) temperature T at different H and (b) 
magnetic field H at different T for Jaa = 1, a = 0.025. 

The temperature dependence of magnetic susceptibility observed in [M] for 
[Ni(H2L^)]4[Cr(CN)6]50H ■ I5H2O compounds resembles the result shown in figure |H 

Notice, that at high temperatures the external field dependence of the magnetic 
susceptibility [figure EJ^b)] exhibits one peak. With decreasing temperature, two 
symmetric peaks begin to arise, which correspond to the formation of the incipient 
(magnetization) plateau at rria = 1/6. With decreasing further temperature the peaks 
become sharper and bigger in their magnitude. 



3.3. Specific heat 

The internal energy u and the specific heat c(T) per cluster site are, respectively, 
determined as 

u = -T^ A (Fcb/^T) (19) 
du T d'^FcB 
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Fgb taken from (fT^ . 




0.000 0.005 0.010 0.015 0.020 



T 

Figure 6. The zero-field specific c(T) heat versus temperature T for Jaa = 1, 
a = 0.025. 

The behavior of the specific heat in equihbrium at the absence of the external 
magnetic field {H = 0) is shown in figure El In this plot one can find the presence 
of second order phase transition: at the same temperature T^, described in last two 
subsections, the specific heat has discontinuity. 
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Figure 7. Specific heat c(T) versus temperature T for Jaa = I, a = 0.025 and 
different H values. 

The temperature dependence of the specific heat at non-zero magnetic field is shown 
in figure [71 Notwithstanding the observation of one peak in the c(T) at if = the 
temperature dependence of specific heat at H exhibits two-peak behavior peculiar 
to one and quasi one dimensional systems as one can find in [551 [56]. A double-peak 
structure in the specific heat manifests the existence of two energy scales in the system 
as a result of two competing orders [57]. Upon increasing the external magnetic field the 
peak moves to higher-temperature region and, at the same time, decreases in amplitude. 
At higher values of H close to the transition from = 1/6 to saturated = 1/2 
state the second broad peak gradually increases. 

In section [SI we further discuss magnetic properties by comparison with thermal 
entanglement. 
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4. Concurrence and thermal entanglement 



The mean- field like treatment of ([H) transforms many-body system to reduced "single" 
cluster study in a self-consistent field where quantum interactions exist. This allows to 
study, in particular, (local) thermal entanglement properties of a-sublattice in terms of 
three-qubit XXX Heisenberg model in effective magnetic field 7a, which carries the main 
properties of the system. Besides, because of the self-consistency and interconnection 
of the fields ja and 7^ the effective 7^ field have an impact on the concurrence, too. We 
study concurrence C{p), to quantify pairwise entanglement [49], defined as 

C(p) = max{Xi -X2-X3- A4, 0}, (21) 

where Xi are the square roots of the eigenvalues of the corresponding operator for the 
density matrix 

p = Puia\^ay)pUa\(^ay) (22) 

in descending order. Since we consider pairwise entanglement, we use reduced density 
matrix pi2 = Trap. Before introducing the calculations and discussion we would like 
to emphasize the fact which was already discussed in section [D the states of two 
neighboring a-type trimers are separable (disentangled). Hence we can calculate the 
concurrence for each of them on cluster level individually in effective magnetic field. In 
our case the density matrix has the following form 

1 ^ 

p = —J2''M-Ek/TMk){^k\, (23) 



k=l 



Ek, \ipk) and Zoa are taken from equations (7), (8) and (9) respectively. The construction 
process of matrix fl22l) does not depend whether 7^ is effective or real magnetic field, 
although the presence of effective field 7^ plays crucial role for the self-consistent 
solution. Here we skip the specific details and provide the result of final calculations of 
the matrix pi2, taking into account that the Hamiltonian "Hco is translationary invariant 
with a symmetry [5*2, "Hcq] = {S^ = Ylk=ii^ z)ki) ■ Hence |i58] : 

/ u \ 



V 



P12 



w 



y 






y 

w 




(24) 



where 



1 



u 



w 



3 
1 

r 
1 

: -e 
3 

1 



27a — 3Aaa 

e 4T 



3(27a+Aaa) 
4T 



4T 



„ 7a „ SAgg 

3e T + 2e 2t 



7g „ 27a+3Aag 

3 + e T -I- 2e 2t 



27g+3Aaa 



y = — e 
^ 3 



1 + e T 



1 + e T 



1 + 2e~ 



-1 + e 2T 



(25) 
(26) 
(27) 
(28) 
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Pi2 in equation is a special case of so called X-state [SSj- The concurrence C(p) of 
such a density matrix has the following form |60j : 

2 

C{p) = -max{\y\ - 0). (29) 

Zj 

Finally, we consider transcendental equations f llSap and fllSfep by taking into account the 
values of variational parameters: \aa = Jaa, Ja = ^Jabrub + H, 'jb = ^Jab'ma + and, 
therefore, one can use these parameters to calculate C(p). First, we study the behavior 
of C(p) at if = 0. The temperature dependence of C(p) is shown in figure [HI 

0.30 
0.25 
0.20 
O 0.15 
0.10 
0.05 
0.00 

0.000 0.005 0.010 0.015 0.020 

T 

Figure 8. Concurrence C{p) versus temperature field for Jaa = 1, a = 0.025 
and if = 

Notice, the " triangle- in-triangle" system can display (bipartite) entanglement 
described by concurrence even in the absence of external magnetic field. It is important 
to mention that this result does not contradict to the well-known fact that there is no 
(bipartite) entanglement, measured by concurrence in isotropic three-qubit Heisenberg 
XXX model in zero magnetic field {H = 0) [H]. Indeed this effect is due to the 
existence of Ising-type interaction replaced by effective field 7^ = 2Jabmb + H acting 
upon a-spins. The latter, in addition to H contains another quantity having meaning 
of magnetic {2Jabmb), which is non-zero at ii = 0. 

Another important observation is that threshold temperature at which entangle- 
ment C(p) disappear is identical to the critical temperature of second order phase 
transition between ordered and disordered phases described earlier in section 13.31 This 
implies that the concurrence vanishes precisely at Tc, the same temperature of specific 
heat discontinuity. This is the consequence of the fact that at Tc the system undergoes 
order- disorder phase transition and the second term in 7^ vanishes, too {nib = 0, when 
H = and T > Tc). This factor implies the strong relationship between magnetic and 
entanglement properties of the system. In figure |9] we present the three dimensional plot 
of the concurrence as a function of the temperature and external magnetic field. We will 
discuss some of these features in behavior of concurrence C{p) for studying magnetic 
and entanglement thermal properties in section |5l 
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5. Common features of magnetic properties and entanglement 
5.1. Finite temperatures 

To the best of our knowledge, the common thermal features in entanglement and 
magnetic properties in frustrated systems have not yet been reported, neither 
theoretically, nor experimentally. In this section we discuss some similarities of magnetic 
statistical properties and quantum entanglement. 

First, we consider for general Jaa and H the susceptibility (fT7|) as a statistical 
characteristic. Figure HU] (a) shows the density distribution of susceptibility reduced 
per one a-site as a function of the coupling constant Jaa and the external field if, at 
a relatively high temperature T = 0.1, which is higher than T^. The white stripes 
on the figure correspond to peaks of the susceptibility. These stripes have a certain 
finite width due to nonzero temperature. For consistency in figure [10] (b) a similar 
plot of concurrence density is shown for the same values of Jaa — H parameters. The 
existence of entanglement in the infinite XXX Heisenberg chains of spins-1/2 and spins- 
1 was pointed in [61] . Here weak probe fields have been aligned along three orthogonal 
directions {x, y and z), supposing that magnetic susceptibility is equal in all these 
directions {xx = Xy = Xz) and using the fact that + + Xz is an entanglement 
witness. In our case there is a magnetic field aligned in the z-direction only, therefore 
Xz = Xa- And we show that the behavior of susceptibility Xz is similar to that of 
bipartite entanglement. Indeed, comparison of figures fTOk and fTOb shows that the 
general behavior of the statistical and entanglement properties, such as susceptibility 
(Xa) and concurrence (C(p)), coincide. Our calculations show that the values of variables 
for the maximum (peak) in magnetic susceptibility correspond to the critical values on 
the Jaa~H diagram at which the quantum coherence disappear and concurrence vanishes 
{C{p) = 0). However, this picture for the Ising- Heisenberg model on the TKL lattice can 
applied only for antiferromagnetic coupling Jaa > 0, while for ferromagnetic coupling 
Jaa < the system always remains disentangled. 
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I 



Figure 10. Density plot for (a) susceptibility Xa and (b) concurrence C{p) 
versus magnetic field H and coupling constant Jaa at a = 0.025 and T = 0.1. 

For further comparison we show for various temperatures in figure [H] (a) and (b) 
the corresponding dependencies of the concurrence and heat capacity on magnetic field 
H. 
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Figure 11. (a) Specific heat c(T) and (b) concurrence C{p) versus external 
magnetic field H for Jaa = 1, a = 0.025. 

In figure [11] (a) at relatively low temperatures the specific heat exhibits six peak 
structure located symmetrically with respect to the magnetic field {H = 0). As 
temperature increases, the middle peaks (on both sides of the H = 0) splits and merge 
with the left and right peaks in the neighbor areas, near H = 0. At higher temperatures, 
the two other peaks on both sides of if = also merge in one. As temperature increases 
there remain only two peaks, i.e., the sharp peak structure gradually disappear. 

At low temperatures close to (T > Tc), the two most distant peaks from the 
H = (on each side of H) in figure [11] (a) are approaching closer to each other, but 
do not merge in one. Meanwhile, the closest ones to the H = peak (on either side 
of H) becomes narrower and approach closer to the origin, H = 0. For T < T^, 
some features are resulting from the effective Ising field: the local minimum of the 
curve c(T) at if = becomes non zero, c{T)\h=o 7^ 0. With further decreasing 
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temperature, the heat capacity in the vicinity of if = displays one narrow peak 
(whereas in a normal three-qubit Heisenberg model there are two symmetrical peaks and 
c(T)\h=o = 0). Nevertheless, the described extreme effects do not affect the behavior of 
the entanglement (concurrence), which is a purely quantum feature: the curve is smooth 
enough for all T > Tc. However, C{p) behaves as a step-like function by approaching 
closer to critical temperature T^. At temperatures below the critical value, a local dip 
minimum is non zero, C{T)\h=o 7^ (see also figure[8]). This dip disappears at T ^ 0.002 
by forming a single flat plateau at T = 0. 

Our effective field results for coupled spins controlled by external magnetic field 
provide understanding of the magnetic ground state properties in " triangles- in-triangles" 
Cu9X2(cpa)6 Kagome series in [62], which are similar to the one dimensional systems. 
This simple approach can also explain several interesting properties, such as double 
peaks of the specific heat, different competing (magnetic) orders, a 1/3 magnetization 
plateau and susceptibility peaks for the pulse field reported for classical and quantum 
Kagome lattice magnets in [63]. In the end of the section, we emphasize that although 
introduced (effective) self-consistent 7a and % fields break the symmetry against H = 0, 
this does not act upon concurrence and specific heat. 

5.2. Zero temperature entanglement and modulated phases 

In this section the magnetization and entanglement properties of a-sublattice are 
considered at zero temperature using variational mean-field approximation. In 
figure [T2r a) a phase diagram of constant magnetization is shown for a-sublattice. This 
diagram differentiates the following phases: Phase I corresponds to the single a-site 
magnetization = 1/6, when spins in a-sublattice are in ffj, configuration; Phase II 
corresponds to H'^ configuration with the single a-site magnetization rUa = —1/6. These 
phases exist only for the antiferromagnetic coupling, Jaa > 0. For the ferromagnetic 
case {Jaa < 0) in III and IV regions we get spin saturation, with maximum rua = 1/2 
{Hi) and minimum ma = —1/2 (ttt) magnetization per atom respectively. 

Phase I contains the two-fold degenerate states {ip^) and {ipe), while Phase Il-the 
two-fold degenerate states \1IJ2) and \ip3)- By constructing the reduced density matrix 
in Phases I and II, one can find these phases in maximum entangled state, C{p) = 1/3. 
Phases III and IV correspond to \ipi) and \ips) states respectively. These phases are 
disentangled, C{p) = 0. In figure fT2](b) the concurrence density distribution is shown 
versus coupling constant {Jaa) and magnetic field {H) at zero temperature. The area of 
non-zero entanglement coincides with the phase I + 11, where \ma\ = 1/6, while the one 
with zero entanglement {C{p) = 0) corresponds to the phase III + IV with \ma\ = 1/2. 

Notice, the plateau behavior in magnetization corresponds to constant entangle- 
ment values in corresponding density plots. The plateau in magnetization at \ma\ = 1/6 
corresponds to maximum entanglement value, C{p) = 1/3, where the saturated phase 
at |ma| = 1/2 is disentangled, C{p) = 0. This descriptive picture is also available at 
the non-zero temperature. At relatively low temperatures the plateau of magnetization 
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at \ma\ = 1/6 and entanglement coincide, except the narrow region in the vicinity of 
H = border. By decreasing the temperature the middle stripe in figure [TOT b) narrows 
and gradually disappears (see figure [T2](b) and also section ISTTl) . This trend becomes 
apparent by comparison of figures [TOT b) and [131 The latter represents the density distri- 
bution of rua magnetization at considerably low temperature T = 0.1. In figure [13] the 
grey areas describe the plateau at \ma\ = 1/6, while black and white regions correspond 
to saturated states, \ma\ = 1/2. White regions in figure [T3] correspond to the plateau 
behavior in the concurrence. As the temperature increases the borders between distinct 
(different) phases are gradually smeared out. Summarizing, the structure of each of the 
Heisenberg trimers has the crucial impact on the phase diagram in figure 12(a): the 
geometrical structure of a-sublattice is responsible for the frustration effects arising in 
the antiferromagnetic Heisenberg model (geometrical frustraion). This leads to above 
mentioned ground states with definite values of concurrence in figure 12(b). The ground 
state concurrence arises on magnetization plateaus at |ma| = 1/6 (see figure 2(b) for 
non-zero temperatures), which is a consequence of strong geometrical frustration of a- 
sublattice. While the octahedron environment (b-sites) are responsible for the effective 
field only. 

(a) (b) 



X 



IV(m„=l/2) / I(m„=l/6) 



III (m„=-l/2) \II (m,=-l/6): 




Figure 12. (a) Phase diagram of a-sublattice for \a\ = 0.025 and (b) density 
plot of concurrence C{p) versus magnetic field H and coupling constant Jaa for 
0.025 at zero temperature. 
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Figure 13. Density plot of magnetization versus magnetic field H and 
coupling constant Jaa for a = 0.025 and T = 0.1. 

6. Conclusion 

In this paper we find strong correlations between magnetic properties and quantum 
entanglement in spin-i Ising-Heisenberg model on triangulated Kagome lattice, 
which has been proposed to understand a frustrated magnetism of the series of 
Cu9X2(cpa)6 • nH20 polymeric coordination compounds. The ratio a = Jab/ Jaa = 
0.025 {Jaa labels intra-trimer Heisenberg, while Jab monomer-trimer Ising interactions) 
is considered, which guaranties experimental realization for suitable theoretical 
treatment. We adopted variational mean- field like treatment (based on Gibbs- 
Bogoliubov inequality) of separate clusters in effective interconnected fields of two types 
(consisting of Heisenberg a trimers and Ising-type b monomers). Each of these fields 
taken separately describes not only corresponding (a- or b- type) spins, but the whole 
system. 

The calculated magnetic and thermodynamic properties of the model display the 
(smooth) second order phase transition from ordered into disordered phase driven by 
temperature and magnetic field. The thermal entanglement properties of a-trimers 
are the central in the aforementioned approximation. Since there are some open 
questions in the definition of multipartite entanglement (for example, see [HUES]) we 
used concurrence as a computable measure of bipartite entanglement for the trimeric 
units in terms of XXX Heisenberg model in self-consistent magnetic field applied 
to a subsystem. Due to the classical character of Ising-type interactions the states 
of two neighboring Heisenberg trimers are separable. Using this fact we studied the 
entanglement of the each "a subdivisions" of the cluster individually in effective Ising- 
type field. Our results show that entanglement does not vanish in zero external field 
as it happens for the mean field treatment of the isotropic XXX Heisenberg model 
on triangulated lattice. The geometrical structure of the lattice is responsible for 
the frustration effects arising in the antiferromagnetic Heisenberg model (geometrical 
frustration). It leads to existence of magnetization plateaus at \ma\ = 1/6 and 
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concurrence at C(p) = 1/3. Besides the effective field in mean-field solution also depends 
on the lattice structure (via the number of 6-sites in the cluster). 

In addition, the entangled-disentangled phases in concurrence and order-disordered 
phases in quantum phase transitions share many common features. The threshold 
temperature for concurrence is identical to critical temperature of second order phase 
transition. Besides the entanglement and thermodynamic properties exhibit also 
common (plateau and peak) behavior in magnetization, susceptibility and concurrence. 
Moreover, we found that for antiferromagnetic interaction the magnetic susceptibility 
peaks coincide with the corresponding phase boundaries at which the entanglement 
vanishes. However, this does not take place for the classical ferromagnetic interaction. 
This fact allows one to notice a quite visible correlation for the boundaries between 
various phases for entanglement, susceptibility and magnetization densities as a 
fingerprints of corresponding quantum phase transitions. Thus, the density diagrams in 
the presence of magnetic field can be considered as a useful tools to detect important 
relationships between entanglement-disentanglement transitions in concurrence and 
corresponding order-disorder quantum phase transitions in quantum magnetism. 
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